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Abstract 

In this work, average Hamiltonian theory is used to study selective excitation in a spin- 1/2 system evolv- 
ing under a series of small flip-angle #— pulses (6 <C 1) that are applied either periodically [which corre- 

o 

^SJ sponds to the DANTE pulse sequence] or aperiodically. First, an average Hamiltonian description of the 

DANTE pulse sequence is developed; such a description is determined to be valid either at or very far 
H- S 

from the DANTE resonance frequencies, which are simply integer multiples of the inverse of the interpulse 

(N 

delay. For aperiodic excitation schemes where the interpulse delays are chosen pseudorandomly, a single 






resonance can be selectively excited if the ^-pulses' phases are modulated in concert with the time delays. 
Such a selective pulse is termed a pseudorandom-DANTE or p-DANTE sequence, and the conditions in 
which an average Hamiltonian description of p-DANTE is found to be similar to that found for the DANTE 
sequence. It is also shown that averaging over different p-DANTE sequences that are selective for the same 
resonance can help reduce excitations at frequencies away from the resonance frequency, thereby improv- 

^ ing the apparent selectivity of the p-DANTE sequences. Finally, experimental demonstrations of p-DANTE 

£T) sequences and comparisons with theory are presented. 
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I. INTRODUCTION 

Of the multitude of radiofrequency (RF) schemes used for exciting and controlling spin dy- 
namics in NMR, most can be placed into one of two categories: aperiodic RF pulse sequences or 
periodic RF pulse sequences. For many aperiodic sequences, the RF phases, amplitudes and pulse 
delays are often chosen randomly or in a pseudorandom manner. Such sequences have been used 
to generate white noise or broadband excitation in NMR noise spectroscopy 1 4 , while sequences 
that generate colored noise have been used in early spin decoupling schemes, such as in noise 
decoupling^. Theoretical models of a spin system's response to pseudorandom pulse sequences 
typically use a Volterra or perturbation series in the randomly applied RF pulses 6 . Since many 
pseudorandom sequences are designed by considering only the first term in the Volterra series, 
pseudorandom sequences are typically low power and result in small, linear spin excitations. 

Unlike aperiodic sequences, periodic RF pulse sequences are commonly used in a variety of 
NMR experiments and are often found to be superior to their pseudorandom counterparts; for ex- 
ample, two periodic sequences, MLEV 7 and WALTZ- 16^, provide better heteronuclear decoupling 
over noise decoupling under most conditions. Many periodic RF pulse sequences are designed us- 
ing average Hamiltonian theory (AHT) 9 , where the necessary RF pulse sequence that generates a 
desired average Hamiltonian H avg over a time r c must be determined (r c is the length of the pulse 
sequence). Repeated application of the pulse sequence introduces frequencies into the dynamics 
that are integer multiples of — , which may result in higher-order contributions to H avg that degrade 
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the sequence's performance. It has been previously noted that random or asynchronous pulse im- 
perfections placed into pulse sequences can often improve their performance 10 . Recently, Uhrig 
dynamical decoupling (UDD) sequences 11 which utilize unequally spaced tx— pulses, were shown 
to be superior in preserving spin coherence to the standard Carr-Purcell-Meiboom-Gill (CPMG) 
sequence 12 , which uses equally spaced n— pulses. 

Selective pulses 13 are one class of pulses that do not fall neatly into either category. The design 
of most commonly used selective pulses, such as the gaussian and the sine pulse shapes, is guided 
by the fact that a spin system's response to an applied pulse as a function of frequency/offset 
is proportional to the Fourier transformation of the applied pulse^. Using linear response to 
design selective pulses has been used to develop colored noise sequences for selective excitation 
in imaging applications^. While sequences designed from the linear response are valid for small 
flip-angles, these pulses fail as the degree of excitation increases. As such, most methods for 



designing selective pulses of arbitrary flip-angles use the linear response pulse shapes as starting 
points in numerical searches 15 . However, pulse shapes generated by these numerical techniques 
often do not lend themselves to an easy physical interpretation behind their selectivity. 

One of the earliest and most easily understood periodic selective pulses that is rigorously valid 
for all flip-angles is the DANTE sequence 17 , which is shown in Fig. [iJA). The DANTE sequence 
consists of a series of N small-tip, broadband 0-pulses that selectively rotate those spins resonating 
at integer multiples of the interpulse delay by 6 = NO about an axis in the transverse plane. The 
DANTE sequence's periodicity is responsible for this frequency response, which can be calculated 
analytically^. To excite a single frequency, however, the periodicity of the DANTE sequence must 
be violated. Breaking DANTE' s symmetry for removing excitation at other frequencies has been 
previously accomplished by modulating the phases^, amplitudes, and delays of the 0— pulses^. 
However, these excitation sculpting modifications of the DANTE sequence are still based on the 
assumptions of linear response. 

In the following paper, we use AHT to provide insight into the selective excitation of a spin- 
1/2 system by a series of periodically and aperiodically small-flip 0-pulses. First, the conditions 
where an AHT description of the DANTE pulse sequence is valid is determined. Next, an AHT 
description for a series of aperiodically spaced and phase-modulated #— pulses is developed. Such 
sequences are referred to as pseudorandom-DANTE or p-DANTE selective pulses [Fig. [jjB)]. 
Finally, experiments performed in acetone and in an acetone/DMSO/water solution are used to 
demonstrate and validate the selectivity of the p-DANTE sequences. 

II. GENERAL THEORY 

Both the DANTE [Fig. [TjA)] and the p-DANTE [Fig. [TjB)] pulse sequences involve the 
application of a series of N small flip-angle #— pulses that selectively rotate spins about an axis 
lying in the transverse plane by an angle = NO. For the DANTE sequence, spins resonating at 
v z = - are selectively rotated by (where n is an integer), whereas for the p-DANTE sequence, 
only those spins resonating at v z = vq are rotated by 0. To understand the selectivity of both 
the DANTE and p-DANTE sequences within the framework of AHT, it is useful to revisit the 
dynamics of a spin- 1/2 system under a non-resonant RF irradiation. The Hamiltonian during the 
application of an RF pulse is given by: ^ = u z I z + ujrf Ux cos(0) + ly sin(0) j, where and 
u RF axe the phase and amplitude of the RF pulse, and u z = 2kv z is the resonance offset that the 
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FIG. 1 . Pulse sequence for (A) DANTEP-^and (B) pseudorandom-DANTE or p-DANTE selective excitation. 
(A) The DANTE sequence consists of a series of N small-flip, #-pulses, equally spaced between periods 
of free evolution of time r. The periodicity of the sequence results in a rotation of 6 = NO about an 
axis perpendicular in the transverse plane for those spins resonating (in the rotating frame) at a frequency 
vz = 7 where n is an integer. (B) The p-DANTE sequence consists of a series of N unequally spaced 
small-flip, 0-pulses where r^ is the time delay between the k th and (k + l) th , and cpk = 27r J2k=i u o T k is the 
phase of the k th pulse with <p% = 0. In the p-DANTE sequence, only those spins resonating (in the rotating 
frame) at a frequency given by vq are rotated by an angle = NO. 

spin experiences in the rotating frame. The propagator for an RF pulse applied for a time T p can 
be written as: 



P^ act (T p ) = exp I -i^Tp J = exp \-i uj z I z + u RF [T x cos(0) + T Y sin(0) 
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Alternatively, the propagator in Eq. ([TJ) can be transformed into an interaction frame defined by 
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ujzlz and is given by: 



PJj xact (T p ) = exp (-iuj z T p I z ) f exp l-i f " dt'^- [l+e i( -" zt ' ~® + I_ e -**>**-*) 

= UfrU"zT p )f exp Li £' dt ^^'A (2) 

where T is the Dyson time-ordering operator, Uf ree (u z T p ) = exp (— iu z T p I z ), and H INT {t'), the 
Hamiltonian in the interaction frame, represents a purely phase-modulated RF pulse. The time- 
dependent propagator in Eq. ([2]) can be approximated by: 

f exp L % - £" dt'H INT (t')\ = exp (-^H p> ^ (3) 

In Eq. Ob, Hp^ is the average Hamiltonian^ and is given by H Pt ^ = Y^=i H p(t> where the first 
two terms in the series are: 

I dt'H mT Jt') 
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P^ xact (T p ) in Eq. (|2|) can be approximated as: 

P?*(T P ) « P^(T P ) » exp (-£Z> Z / Z ) exp (-^ (s£J + B*J)) (5) 

For tu RF T p < f , ||P^ xact (T p ) - p0(T p )|| < 10 -3 for all w z , where ||A|| = y / Tr[AUf represents 
the Frobenius matrix norm (if A represents the difference of two unitary matrices, then the max- 
imum value of \\A\\ is 2y/n where n is the matrix dimension). Since the DANTE and p-DANTE 
sequences both consist of a series of small flip-angle #— pulses with 6 < ^, the approximation 
P| xact (T p ) « P^iTp) in Eq. (pb will be used in the rest of this paper. 

For future comparison of the propagator in Eq. ([2]) with the propagator for a spin- 1/2 evolving un- 
der either the DANTE or the p-DANTE sequences in Fig. [T] it is useful to consider an alternative 
description of P^ xact in the interaction frame by dividing Texp ( — | / p ifi NT (t / )dt'J in Eq. Q2b 



into iV ^> 1 smaller propagators, which is illustrated in Figure [2?A). In this case, P^ xact (T p ) can 
be rewritten as: 

,HlNT(t')\ 
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« U^(u z T p )T f[ exp (-t^psinc (^) (T x cos (</>*) - I y sin (</>*))) (6) 

where 0* = "^ - p + ^ — 0. With respect to Eq. (6 ) and Fig. 2[A), the total propagator for an 
RF pulse of strength curf applied off-resonance by coz for a time T p is equivalent to the application 
of N continuous, small-flip 6 = ^sinc (^r), phase-modulated RF pulses where the phase of 
the j th pulse is —(j)% followed by a rotation about the z— axis by u>zT p . 

A. DANTE Pulse Sequence 

The DANTE sequence 17 consists of a series of N, equally spaced small-tip, ^-pulses of constant 
phase and duration t p where ujRFt p = 9 <^ 1 [Figure [1? A)]. The full propagator for the DANTE 
pulse sequence can be written as: 

t4xact(Ttot) = Po(tp) (Ufve S (UzT)P (t p )) = (Pofe)^free(w 2 r)) P (t p ) 

N-l 

= UiUuj z ((N - l)r« + t p ))f [] P-k. Z T t (t P ) (7) 

fc=0 

where r is the time delay between pulses, r t = r + t p and 

P-kuzn(tp) = UL^zkr^Poit^Uf^iuzkn) = exp (-iH k t p ) (8) 
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(9) 

As has been previously notecP^H the propagator for the DANTE sequence in Eq. ([7j) is the 
same as the propagator for a continuous series of N, phase modulated small-flip pulses, where the 
phase modulation depends upon the spin's chemical shift, uz, followed by a rotation about the 
z— axis by u>zT m . This is illustrated in Figure|2|B). If mod [u>zTt, 2ir] ps 0, then all N pulses are 
effectively applied along the same direction since 4> k « <pj for all k and j. The small rotations 



are therefore additive and lead to an overall rotation of 9 ~ NO about an axis in the transverse 
plane is generated. When mod [uj z r t , 2n] 7^ 0, the pulses are effectively applied about different 
directions (0 fc 7^ <pj for k 7^ j in general) thereby reducing the overall spin rotation. Comparing 
Fig. gB) and Eqs. and ^ with Fig. [2f A) and Eq. ((6]), the propagator for the DANTE sequence 
is similar to that of an off resonant, RF pulse of duration Nt p = T p followed by a rotation about 
the z-axis. That is, 

U e xact(Tm) ~ Uhe S (uzTtot - uj' z Nt p ) exp -iT p (u' z I z + u' RF (Jx cos(0') + ly sin(0'))) 

(10) 

where u' z = |^mod [cuzTt, 27r], u' RF = : — \, 2 t ( corf, and 0' = -f [uj' z — u z ). As mentioned 

above, when mod [cu z r t ,27r] = 0, then co' z = 0, and the effective pulse is applied on resonance 
and rotates the spin by 6. When mod [u Z T t , 2n] 7^ 0, then uj' z can be quite large since j- 3> 1 for 
short pulses (t p <C 1). In this case, the pulse appears to be applied very far off resonance when 
u' z > u' RF . 

In order to make the above arguments more quantitative, AHT can be used to rewrite U(T tot ) 
in Eq. (J7J as U(T tot ) « U Am (T tot ) where U Am (T tot ) = U fme {u z T tot )ey^(-iNt p H av ^, where 
the first two terms in the average Hamiltonian, H avg = Y^Li H avg , are given by [setting a = 
cuisine ( ^ ) and b = %£ (1 - sine (u z t p ))]: 
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where f T (u z , r t , £ p , N) = T x cos ( i± ^ k ) - T Y sin (^f^) . Although the form of H avg in Eq. ( 1 1 ) 



is somewhat complicated, the physical picture behind H avg in Eq. (11) can be seen in Fig. \2 B) 
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When u Z T t = 2ixn for integer n, <pj = <pk and therefore Hk, Hj = for all k 7^ j, so H 
H a J exactly. In this case, the propagator represents a rotation about an axis in the transverse 
plane of phase (—l) n co z t p /2 by a total angle of 6 = NuiRFt p smc (u) Z t p /2) t=a NujRFt p = N9 
for uj z t p <C 1. For uo Z T t 7^ 27m for integer n, the various H k are pointing in different directions, 
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FIG. 2. The connection between the DANTE pulse sequence and the application of an off-resonant RF 
pulse. (A) An RF pulse applied off resonance by ojz and with strength ojrf for a time T p is equivalent to 
applying N S> 1 phase-modulated, small flip-angle = Wiyr-^sinc ( ^^n" ) P u l ses followed by a rotation 
about the £-axis by an angle ojzT p . The phase of the j th pulse is given by —<fjj = (j — l)ojz-fi +wz^ — 
(in Fig. [2]^A), = 0). In (B), a DANTE pulse sequence is equivalent to a series of N phase-modulated 
small-flip pulses followed by a rotation about the z— axis by an angle Nujzt [Eq. Q], The phase of the 
k th pulse is </>fc = ujz \(k — l)(rt) + -f ) where T t = r + £ p , t p is the length of the small-flip pulses, and 

T tot = {N- i)r t + t p . 



so that the average transverse field in H is lessened. Furthermore, since the effective rotation 
directions no longer commute with one another, i.e., [H k , Hj] ^ for k ^ j, there is a contribution 
to H avg at second-order, H avg , of an effective field along the z direction. When \H avg \ 3> \H \, 
the effective field lies mostly about the z— axis and the spins are minimally excited (this argument 
is similar to the concept of second- averaging^). 

In order to see under what conditions AHT can be used in calculating the DANTE pulse se- 
quence, Figure |3J shows 1 1 C/ ex act(7tot) - ^AHr(^tot)|| using # avg « H^ g + H^ g , as a function 
of N and luz/urf- Two calculations are shown in Fig. [3j one for a total pulse rotation for 
an on-resonant RF pulse of = n/2 [Fig. [3^ A)] and one for = 2% [Fig. J3^B)]. For both 
calculations, f- = 1000. As mentioned above, U A uT(T to i) ~ C4xact(^tot) when 2nn = tuz T t = 

^~n \t~ ~*~ ■*■)' wnere n i s an integer. Therefore, at a resonance condition 
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for n^O, there exists a linear relationship between N and -^- that is given by 
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Urf 2vrn \t p 

where INT[x] gives the nearest integer near x. From Fig. [3j for 9 < ~ [N > 30 in Fig. [3VA) 
and N > 120 in Fig. 3 13)], UAm(T t ot) is a good approximation to U exact (T tot ) for all ^- except 



near the resonance conditions in Eq. (12). The deviations of t/AHT^ibt) from t/ exac t(T tot ) occur 
when u)z is slightly away from the resonance condition, u z = — for integer n. From numerical 
calculations, the range of frequencies in which £/AHr(^Tot) is a good approximation to U exact (T t0 ^) is 
found to be approximately given by 5ui z 3> J^ or 5u z <C J^, where 5u z = min \u z — ^p| 
is the smallest frequency difference between co z and the nearest integer multiple of the DANTE 
resonance frequency, — . Although the range of u z where Uxm:{T lo ^) is a good approximation 
increases with decreasing 9 (increasing N), the maximum of \\UAm{T t ot) — ^exact(^toi)ll mainly 
depends upon the overall rotation angle, 0. For 6 < ^, max \\U\wr{Ttot) — U exact (T t 



-tot 



< 



0.1 for all iV and co z . Physically, this can be understood as follows: for 5u z t p ^$> J^, the 
effective phases of the pulses [see Fig. [2p3)] are modulated faster than the effective tip of the 
pulse, URFtp = 9, so that co z ^> corf- In this case, AHT works well, just as in using the typical 
rotating wave approximation. For lo z <ti J^, the phases of the pulses in Fig. 2TB) are relatively 
unchanged during the course of the sequence; in this case, the various Hamiltonians, Hk in Eq. (|9|) 
commute with one another, so H avg w H* vg and U^ act (T tot ) w U Am (T tot ). 



B. pseudorandom DANTE (p-DANTE) 

In the DANTE sequence, a natural frequency of — is introduced into the dynamics due to the 
periodicity of the pulse sequence; this leads to efficient excitation at frequencies z/° ANTE = — 
for integer n. However, suppose that one was interested in using a DANTE-like sequence to 
efficiently excite only one particular frequency, say at Uq ANTE = Hz. One way to accomplish this 
using a DANTE sequence would be to make r t small enough such that all v °^ TB lie outside the 
relevant spectral width. If the spectral width for the system of interest is large, however, this would 
necessitate using small r t , where the smallest time possible r t is r t ps t p (i.e., when r = 0). The 
selectivity or width of the excitation spectrum about z/° ANTE is approximately given by j^- . For 
Tt ~ t p , this means the selectivity is roughly proportional to ^- = 2ia ^ F . In this limit, the effect 
of the DANTE sequence is similar to evolution under continuous RF irradiation, which leads to a 
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FIG. 3. The Frobenius norm of ||t4xact(^tot) — ^AHT(^tot)|| as a function of the number of small tip 8- 
pulses, N, with 9 = ^, and ■^ z ~- In Fig. 3 two calculations are shown for total maximum flip angle of 
(A) = | and (B) = 2ir. C/ exact (T tot ) [Eq. o7|] is the exact propagator for the DANTE sequence, and 
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^AHr(Ttot) is the propagator calculated using the average Hamiltonian, H avg « H y ^ g + H ai , 



in Eq. (JTTj). 



In both calculations, ~ = 1000. The greatest deviations between f7exactC?tot) an d ^AHr(Tt t) occur about 



the resonance conditions in Eq. (12), and the maxima in ||t/ e xact(^tot) — ^AHTC?tot)|| are approximately 
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described by two parallel lines given by iV = INT , ,__ < , 
agreement between UASi^T\at) an d C4xact(?tot) improves with decreasing 0. 



for integer n / 0. The 



very broad excitation profile unless v RF is weak or ^ 2n. Under these conditions, the DANTE 
sequence would be equivalent to applying a long, low-amplitude RF pulse. 

An alternative way to excite only a single resonance using a DANTE-like sequence would 
be to violate the periodicity of the DANTE sequence. This could be accomplished in a variety 
of ways, such as using aperiodic delays, modulating the pulse amplitudes and delays, etc. One 
example of such an aperiodic DANTE sequence is illustrated in Figure [jjB) where N small-flip 
^-pulses are applied with a delay of r^ between the separation between the k th and [k + l) th pulse, 
where in general, Tu ^ Tj. For such pulse sequences to selectively and efficiently excite spins at 
a single resonance frequency u , the phases of the pulses, <p k , must be modulated. In this case, 
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the phase of the (k + l) th pulse is given by fc+1 = — Ikvq (kt p + HjLi T j) = ~ 27rz/ T fc with 
T fc = kt p + Sj=i r j> and 4>\ = and T = 0. Such a set of aperiodic DANTE sequences are 
referred to as pseudorandom-DANTE or p-DANTE sequences. 

As in the DANTE case, the propagator for the p-DANTE sequence can be written as U(T tot ) = 
Ufree(uz(T t ot)) exp (-j.Nt p H avg j , where T tot = T N _ X + t p , and H avg is the average Hamiltonian 
for the p-DANTE sequence, with the first two terms given by: 

N / N 

< = £ E ( v^*^) + /- e -( A - Tfc+ ^)) + & d3) 

avg 2N 

+ ^^£sin(A W (T fc -T J )) < 14 > 
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where Aw = 2n [yz — ^o), and a and 6 were previously defined before Eq. (11). Unlike in the 
DANTE case, the average delay between pulses fluctuates within the p-DANTE sequence, i.e., 
^ ^ -4- for k 7^ j. For vz ^ v®, a p-DANTE sequence is effectively equivalent to applying an 
RF field with a fluctuating offset [Fig. [2^A)], the result of which is a seemingly random excitation 
profile for those spins with Acoz ^ 0. Spins with Aw z = are rotated by N9 = 6. This 
is illustrated in Figures ^A) and[4pB), which show the excitation and z— magnetization profiles 
under a p-DANTE sequence respectively. In Figs. |4jA) and[4jB), the various N — 1 Tfc's were 
chosen randomly but were scaled to ensure that Ylk=i T k — 46.4 ms, and one hundred different 
sets of randomly generated p-DANTE sequences were generated. Consider the excitation and 
z— magnetization profile for a single p-DANTE sequence [red curve (N avg = 1)] shown in Figs. 
41 A) and 413). A maximum rotation by 9 = | occurs at 5v = 0, where — (Iy) = 1 and (l z ) = 0. 
Away from bv = 0, the excitation and z— magnetization profiles are quite noisy, but | (Iy) \ < 1 
and (lz) > 0. Note that the z— magnetization profile is less noisy, since rotations away from the z 
direction go as cos(6>) « 1- y whereas excitations go as sin(0) ~ 9 for 9 <C 1. 

If the excitation and z— magnetization profiles are averaged over different p-DANTE sequences 
that possess the same total pulse length and are all selective for z/ , then the fluctuations in both 
the excitation [Fig. [4jT3)] and z— magnetization [Fig. |4]T3)] profiles for 5v ^ decrease rela- 
tive to the excitation at 8v = roughly as J- — . However, even when N avg ^$> 1, there still 

\J Navg 

exists a "baseline" excitation at 5v ^ which is nonzero (averaging simply decreases the fluc- 
tuations about the baseline excitation). The average "baseline" excitation is approximately given 

11 



by \(Iy)\ ~ ^2 anc ^ ( z ^ = ^ ~ 2 • Thus decreasing 9 and N will decrease the "baseline" 
excitation as 9 oc i. Similar schemes averaging over random sequences have been previously 
used for stochastic dipolar recouplingj 25 l 26 [ 

Besides randomly chosen delays, averaging over different sets of delays that are periodically 
modulated can also lead to selective excitation. Consider a series of delays where the k th delay is 
given by Tk = t + 5t cos ( ^) , where / is a real number, and r > 5t so that Tk > 0. For such a 
sequence to selectively excite spins resonating at u , the phase of the k th pulse must be given by 
(j) k = —2iris T k _i with T = and : 

Tk = }^ T k = kr - — I 1 - esc I - I sin I — - — vr I I (15) 

where the total time of the sequence is given by T tot = Nt p + 7/y-i = Nt p + (N — l)r — ^ + 

f esc (J) sin (^ii). 

Using the values of T fc in Eq. ( 15 ), H av in Eq. ( 13 ) can be evaluated and is given by: 



avg 



*%. - 4 E ( he< A " n+ ^ + 7 _ e -<(^+^) ] w . 



fc=i 






n— — 00 



- \J / ' sine I g — 



where Aw = u;.z — u; , J„ is a bessel function of order n, and Xn = Aw- — -^ — - + ^j 11 + ^f^- 



From Eq. (|16|), if is maximal at the conditions 2it/\vt + "^f- = 2miT or at Av = — — -?- where 



. ^ .uu»iu™ w «v wuwuwuo «-„ «^, I y, ~w„„^.«,— T ^ 



m and n are integers. These define the resonance conditions for this type of p-DANTE sequence. 
However, H is scaled by J n ((m — j) ^ esc (?)), which is greatest when m = n = 0. For 
m^O and n^O, this scaling is less than one, which results in a smaller total rotation. 

Figure [4jQ and (D) show the numerically averaged ^-magnetization and excitation profiles 
respectively, averaged for up to one hundred different p-DANTE sequences using periodically 
modulated delays. In the simulations, N = 30 and ^ = -4=. For the p th p-DANTE sequence, f p 
was set to be equal to the inverse of the square root of the p th prime number, i.e., /1 = l/\/2, f 2 = 
1/V3, /100 — 1/V541. In order to better compare these results to the results for the p-DANTE 
sequences using random delays [Figs. WlA) andWfB)], r v for the p th experiment was chosen so 
that ^ Z)?=i tJ = 1.6 ms. First consider the A^ = 1 case (red curve) in which / = -4= and 
i = 625.13 Hz. Unlike the case of using random delays [red curves in Fig. RlA) andWTB)] where 
the resulting excitations appear randomly distributed throughout the spectral range, the excitation 
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FIG. 4. Theoretically calculated "averaged" excitation [(B) and (D)] and z— magnetization profiles [(A) and 
(C)] for a series of p-DANTE sequences [Fig. [TlB)]. In all p-DANTE sequences, N = 30, 9 = gg, and 
4>k = for all k so that only spins resonating at 5v = Hz are rotated about the x-axis by G = f - In (A) and 
(B), one hundred different p-DANTE sequences with randomly chosen delays were generated such that the 
average pulse delay, ^g Ef=i T fc = 1-6 ms, was the same for all p-DANTE sequences. The excitation and 
z-magnetization profiles were averaged over N avg p-DANTE sequences with N avg = 1 (red curve), N avg = 
25 (blue curve), and N avg = 100 (green curve(. In (C) and (D), one hundred different p-DANTE sequences 
were chosen such that the k th delay for the p th experiment was t% = t p 1 + ^ cos ( ^ J where ^ = 4= 
and f p was inversely proportional to the square root of the p th prime number (e.g., f\ = \j\[2, J2 = l/\/3, 
/ioo = 1/V541). t p was chosen to ensure that the averaged delay was ^ Sl=i T k = 1-6 ms ^ or a ^ P- The 
averaged excitation and 2-magnetization profiles for N avg = 1 (red curve with f\ = l/\/2), N, 



"avg 



25 



(blue curve, averaging from f\ = 1/v 2 to /2s = 1/V97), and N avg = 100 (green curve, averaging from 
/1 = l/\/2 to /ioo = l/sgrt541) are shown. As N avg increases, the "fluctuations" in both (Iz) and —(Iy) 
decrease, and the results are similar for p-DANTE sequences using randomly chosen delays [(A) and (B)] 
and those using periodically modulated delays [(C) and (D)]. 
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profile using periodically modulated delays occur at discrete 8v given by the resonance condition 



8v = 625.13 (m — \/2nj Hz [Eq. (16)]. Note that while the resonance at 8v = (m = and 
n = 0) is maximally excited ((Iz) = and —(iy) = 1). the degree of excitation at other resonance 
conditions is less. In particular, the resonances at 8v = ±625.13 Hz (m = ±1, n = 0) are not 
observed in the calculated profile, since at these conditions, H avg is scaled by Jo( _ 2.3046) = 
0.053, whereas the resonances at 5u = ±366.2 Hz (n = =Fl and m = ±2) are clearly observed 
(| Ji (1.35) | = 0.5325). As was the case for p-DANTE sequences using randomly chosen delays, 
averaging over different sets of periodically modulated p-DANTE sequences reduces the excitation 
for all resonances except at 8v = 0, which is a common resonance for all p-DANTE sequences. 
From Fig. |4} the ^-magnetization and excitation profiles using periodically modulated delays 
[Figs. |4jC) and|4fD)] become similar to those using the randomly chosen delays [Figs. |4jA) and 
WD)] as N avg increases. 

Finally, it should be noted that the conditions under which the average Hamiltonian in Eq. 
( [To] ) provides a valid description ofthe p-DANTE sequence are approximately the same as those 
found for the DANTE sequence [Fig. |3j. Figure [5] shows the difference in the excitation and 
^-magnetization profiles calculated using either the exact propagator or the propagator calculated 



using the average Hamiltonian up to second-order [Eq. 14 and Eq. (16)] for the p-DANTE se- 
quences used in Figs. |4jC) and|4^D). AHT works well for all 8v away from resonance conditions 
which is evident from Fig. 15] for the N avg = 1 curve. The magnitude of the error in this approx- 
imation is the same as that found for a DANTE sequence with = n/2. The agreement of the 
AHT calculations with the exact calculations appears to improve upon averaging over different 
p-DANTE sequences, except for the 8v — resonance. This is due to the fact that only the 8v = 
resonance is the same for all p-DANTE sequence used in Fig. [5} 

HI. EXPERIMENTAL 

All experiments were performed on a 300 MHz Avance Bruker spectrometer (static magnetic 
field of 7 T and an operating frequency for l H of 300. 1 3 MHz), using a 5-mm Bruker BBO probe. 
A 2M solution of acetone in acetone— <i 6 was used to experimentally determine the excitation and 
^-magnetization profiles as a function of frequency offset from the acetone resonance for both the 
DANTE and two pseudorandom pulse sequences. The carrier frequency was incremented between 
—580 Hz below to 580 Hz above the acetone resonance in intervals of 10 Hz in order to exper- 
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imentally determine the excitation and ^-magnetization profiles (a total of 1161 measurements), 
and the integral of the acetone peak was measured. In order to measure the ^-magnetization, a 
| pulse (Rabi frequency of 21.4 kHz) was applied after the DANTE and p-DANTE sequences, 
which was phase cycled in concert with the receiver phase so that only the z— magnetization prior 
to the last | pulse was measured. A delay of 40 seconds was used between scans in all experiments 
in order to ensure that the system had relaxed back to equilibrium which was necessary to avoid 
any distortions in the observed profiles. 

In order to demonstrate the improved selectivity in the excitation and z— magnetization pro- 
files by signal averaging over different p-DANTE sequences (as shown in Fig. [4]), experiments 
using different p-DANTE sequences were performed on a solution of acetone, dimethyl sulfoxide 
(DMSO),and water diluted in D 2 0, such that ^ggJ = 0.822 and Jjggj = 1.4355. All chemicals 
were obtained from Sigma- Aldrich. 

IV. RESULTS AND DISCUSSION 

The experimentally determined excitation and z— magnetization profiles under the DANTE 
and two different p-DANTE sequences obtained using a 2M acetone solution in acetone-rf 6 are 
shown in Figure [6j where the blue and red curves correspond to the theoretical and experimentally 
observed profiles respectively. In these experiments, iV = 30, 6 = ^, and t p = 720 ns were 
used with a maximum rotation of = NO = |. For the DANTE sequence, r = 2 ms. Over 
the spectral range shown in Fig. pi A) andpTD), excitations at frequencies 5v = ^ = ±500 Hz 
and at bv = Hz were observed. For the p-DANTE sequences, the k th delay was given by either 



[Figs. 



6;C)and 



6F)] r k = 2.063 [1 - cos (^ I JJms=4.126sin^ (^f^yj (which is a similar set 
of delays used in the UDD sequences™) or [Figs. Rb) andlefc)] r k = 2.096 Pi + I ms ( ™«\ 



3 COS V 23 j 



In both cases, the average delay between pulses was equal to ^ Y^=\ r k = 2 ms in order to allow 
for better comparison with the DANTE sequence used in Figs. [6jA) andj6jD). Both p-DANTE 
sequences generated a maximum excitation at 8v = Hz, and smaller excitations for 5v ^ 
Hz were also observed, as expected. Note that for the UDD-like p-DANTE sequence [Fig. [6fC) 
and[6]T)], the excitation and 2-magnetization profiles look similar to that of a p-DANTE sequence 
using randomly chosen delays [the red curve in Figs. |4jA) and[4pB)], whereas excitations using 
the other p-DANTE sequence [Figs. [6jB) and[6jE)] appear to be concentrated within a smaller 
frequency range. 
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In order to examine the effects of averaging over different p-DANTE sequences, experiments 
were performed on a DMSO-acetone-water solution in D 2 0. The spectrum of the solution after 
a simple |-acquire sequence is shown in Fig. P7VA), where the RF was applied on resonance with 
respect to the water resonance [<5z/ acetoneiWater = —768 Hz and 5z/ DMSO , W ater = —620.2 Hz]. The 
experimental excitation and z— magnetization weighted spectra after the application of a DANTE 

sequence with N = 30, 9 = Jg (t p = 630 ns), and r ps -^ — r = 1.6 ms are shown in Figure 

|7jB) and Figure [T^C) respectively. With this choice of r, the DANTE sequence efficiently excites 
both the water (Su = 7) and DMSO (Su = — M resonances [Fig. 7 IB)] and leaves the acetone 
magnetization mostly about the z— axis [Fig. |7fC)]. 

The averaged excitation [Fig. [8^B)] and z— magnetization weighted spectra [Fig. [8|A)] for 
the DMSO/acetone/water solution was obtained using up to one hundred different p-DANTE 
sequences, and the results are shown for N avg = 1 (red curve), N avg = 25 (blue curve), and 
N aV g = 100 (green curve) in Figure [8] The p-DANTE sequences used in Fig. [8] were the same 
used in the theoretical calculations shown in Figs. (4j^C) and|4[D), where the k th delay used in p th 
experiment was given by: 

p 46.77ms / 1 /2ttA;\\ 

which ensured that the average delay, ■— Y^Li T l ~ w, — " 1 = 1-6 ms, was the same as the 

delay used in the DANTE sequence shown in Figs. [T^B) and[7]^C). As in Figs. [4jQ and[4jD), 
the water (Su = Hz) was maximally excited whereas the averaged excitation at the acetone and 
DMSO resonances decreased upon averaging over different p-DANTE sequences. Similarly, the 
^-magnetization weighted spectra indicated that the acetone and DMSO magnetization remained 
mostly about the z-axis after application of the p-DANTE sequence, whereas there was little z- 
magnetization at the water resonance. 

V. CONCLUSIONS 

In this work, average Hamiltonian theory (AHT) was used to calculate the effective propaga- 
tors for the both the DANTE [Fig. [TJA)] and pseudorandom-DANTE or p-DANTE [Fig. [TjB)] 
sequences. It was found that an AHT description the DANTE sequence is valid when 9 = ~ and 
for total pulse flip-angles of 6 = iV^ < ^ over all frequencies [Fig. |3l. The validity of the 
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AHT description was also found to depend upon the spin's resonance frequency, v, and an AHT 
description of DANTE works well for frequencies in the range 8v ^$> -^^ and bv <C ^q where 
5u = min 



v 



is the smallest frequency difference between v and the nearest resonance of 
the DANTE sequence, — where n is an integer. Understanding the limitations of an AHT descrip- 
tion for DANTE enabled us to develop an AHT description of the p-DANTE sequences [Fig. [TpB)] 
where the delays and phases of the pulses are modulated in concert throughout the sequence. This 
modulation of delays and phases breaks the periodicity of the DANTE sequence and enables the 
p-DANTE sequence to excite spins at a single frequency, v . The ability to use an AHT descrip- 
tion for the DANTE and p-DANTE sequences might also provide additional insights into other 
selective pulse sequences, since any shaped pulses can be cast into a DANTE-like description^^!]. 
While the excitation and z— magnetization profiles for a single p-DANTE sequence are not partic- 
ularly clean, i.e., small excitations exist at many frequencies, averaging over different p-DANTE 
sequences helps to "clean-up" the excitation profiles so that only a baseline excitation exists ev- 
erywhere except at u , which is excited. Experimental demonstrations [Fig. [6] and Fig. [8j of the 
p-DANTE sequences were found to be in good agreement with theoretical predictions. 

For future work, determining the existence of an optimal set of p-DANTE sequences that gen- 
erate the "cleanest" excitation profiles using the smallest number of p-DANTE sequences will be 
investigated. Since the frequency selection in p-DANTE sequences is determined by correlating 
the pulse phases with the delays, the p-DANTE sequences could also be incorporated into ultra- 
fast NMRpS^H techniques to selectively excite certain resonances in different parts of the sample 
volume. Furthermore, extending the AHT results obtained in this paper to coupled spin systems 
is currently underway, whereby a DANTE-like or p-DANTE-like sequences can be used to selec- 
tively excite a particular multiple-quantum spin transition. The conditions under which such an 
AHT description can be applied in these systems are approximately the same as those found in this 
paper, since any subspace of two transitions can be described 29 as an effective spin- 1/2. Coupling 
these techniques with ultra- fast NMR should enable the quick determination of all spin transitions 
in a given molecular system. 
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FIG. 5. Difference in the calculated (A) z-magnetization and (B) excitation profiles for the p-DANTE 



sequences used in Fig. |4|C) and WD) between the 



UAKr(Ttot) 



iT^ z T z - iNt ^l+<) 



profiles calculated using the propagator from AHT, 



|l4 



[Eq. (16 1 and Eq. dl4i] and the profiles calculated using 



exact propagator, U exact (T tot ). The agreement is relatively good over a wide frequency range except near 



the 5v = Hz resonance condition, 5v 



39 
-50r t 



±13 Hz. For all \Su\ » 



5^ and \6u\ « J^ 



the agreement between AHT and the exact calculation improves with averaging over different p-DANTE 
sequences. 
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FIG. 6. The experimental (red) and theoretical (blue) excitation [(A)-(C)] and z-magnetization [(D)-(F)] 
profiles under application of the DANTE [(A),(D)] and p-DANTE sequences [(B),(E),(C),(F)] as a function 
of the applied RF's offset away from the acetone resonance, 5v. In all experiments, N = 30 and 8 = ^. 
The profiles were generated using a 2M acetone solution in acetone-dg by changing the RF carrier frequency 
from —580 Hz to 580 Hz in 10 Hz increments, and the resulting acetone resonance was integrated. For the 
DANTE sequence, r = 2 ms; with these parameters, the acetone resonance was maximally excited at 



5u 



±500 Hz and at 5u = Hz over the spectral range [-580 Hz, 580 Hz]. For the p-DANTE 



sequences, the k th delay was either given by r^ = 2.096 1 + \ cos ( %_ 



the UDEffH-like delay r k = 2.063 



1 — cos 



kir 
29+1 



ms=4.126sin 2 



ms 



[Figs.0B)and[6|E)]orby 
ji^jy [Figs, go and^F)]. In both 



cases, the average delay between pulses was ^ 2l=i T k = 2 ms in order to enable comparison with the 
results from the DANTE sequence in Figs. (6jA) and[6jD). For both p-DANTE sequences, there is minimal 
excitation at bv = ±500 Hz, although the UDD-like sequence [Figs. (6jC) and|6jF)] appears to generate 
smaller excitations over a wider frequency range than the other p-DANTE sequence [Figs. |6jB) and(6jE)]. 
In all cases, there is good agreement between theory (blue) and experiment (red). 
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FIG. 7. (A) The spectrum after a ^-acquire experiment for a DMSO, acetone, and water solution. The 
spectrum is centered on the water resonance (5u = 0). (B) The spectrum after application of a DANTE 

pulse sequence with TV = 30, 9 = ^, © = NO = |, and r = ^ r w 1.6 ms. In this case, both 

the water and DMSO resonances are excited whereas very little excitation occurs at the acetone resonance 
(<5^acetone,water = —768 Hz). (C) The z— magnetization weighted spectrum after application of the DANTE 
pulse sequence. As expected from (B), there is substantial z-magnetization for the acetone resonance and 
little z-magnetization for both the water and DMSO resonances after application of the DANTE sequence. 
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FIG. 8. Experimental stacked plots of the (A) z-magnetization weighted spectra and the (B) excitation 
spectra after averaging over N avg = 1 (red), N avg = 25 (blue), and N avg = 100 (green) p-DANTE 
sequences [same sequences used in Figs. |4jC) and[4jD)] applied to the acetone, DMSO, and water solution 
used in Fig. [7] From Figs. |4jC) and |4^D), only the water resonance at 8v = Hz should be efficiently 
excited. From Fig. (6lB) the water resonance (5u = 0) is efficiently excited, and the amount of excitation at 
the acetone and DMSO resonances decreases with averaging over more p-DANTE sequences. In[8]A), the 
z— magnetization weighted spectra are shown, illustrating that both the acetone and DMSO magnetization 
lie mostly along the z-direction after application of the p-DANTE sequences. 
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